Bose-Einstein condensation at constant temperature 
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We present a novel experimental approach to Bose-Einstein condensation by increasing the par- 
ticle number of the system at almost constant temperature. In particular the emergence of a new 
condensate is observed in multi-component F = l spinor condensates of *^Rb. Furthermore we de- 
velop a simple rate-equation model for multi-component BEC thermodynamics at finite temperature 
which well reproduces the measured effects. 
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The experimental realization of Bose-Einstein conden- 
sates (BEC) in dilute atomic gases Q, H, |1| and the 
breathtaking emergence of fascinating physics of cold 
quantum gases in an increasing number of experiments 
have had formative influence on the common model usu- 
ally used for the description of Bose-Einstein condensa- 
tion (see e.g. and references therein). This model is 
based on a system of constant particle number whose 
temperature T is reduced. The popularity of this ap- 
proach arises from the fact that all experiments so far 
make use of evaporative cooling techniques which reduce 
the temperature of the sample (at the expense of particle 
losses). This path to quantum degeneracy is illustrated 
in the phase diagram of Fig. ^ Starting with a certain 
particle number N, the temperature T of the system is re- 
duced below the critical temperature Tc{N) which leads 
to an accumulation of particles in the condensate frac- 
tion Nq/N. Detailed experimental studies 0,0,01 have 
compared this quantity with theoretical descriptions. 

In this paper we present a completely different and 
new experimental realization of Bose-Einstein condensa- 
tion by increasing the particle number of a system at 
almost constant temperature. The corresponding path 
is also marked in Fig. ^ and leads to BEC almost or- 
thogonally to the common route discussed above. We 
start with = and add more and more particles at 
nearly constant temperature T until the critical particle 
number Nc(T) is reached, i.e. the population of the ther- 
mal cloud saturates and all further added particles fill 
up the condensate fraction. Worth mentioning this ap- 
proach corresponds to the original idea used by Einstein 
|8j and theoretical descriptions over decades to discuss 
BEC. Furthermore first attempts to achieve quantum de- 
generacy in spin-polarized hydrogen ||9||, 110| were based on 
increasing density by adding particles and by compres- 
sion at liquid helium temperatures. 

The new thermodynamical approach to BEC discussed 
in this letter is realized in multi-component BEC which 
provide multiple internal quantum states of the involved 
atoms. We want to emphasize that these systems open 
up a rich variety of new thermodynamical aspects as the 
involved finite temperature dynamics is extended to more 
components which are additionally coupled and influence 




FIG. 1: Phase diagram of Bose-Einstein condensation for a 
typical *^Rb experiment. The condensate fraction (if > 0) is 
plotted as No/N = 1 - g3{l){kB / (huj)fT'-^ /N . The usual re- 
alization of BEC is done by decreasing T at (almost) constant 
particle number A'^. In this paper condensation by increasing 
particle number starting with = at (nearly) constant 
temperature is discussed. 



each other. The thermodynamical description has to take 
into account all interactions between multiple condensate 
components and just as many thermal clouds (we use this 
term instead of 'normal components'). In this context re- 
cent experiments have observed 'decoherence-driven cool- 
ing' and melting of new condensate components ■ 

The system considered here is based on a = 1 
spinor condensate of ^^Rb with three internal states 
mp = — 1,0,+1. The main idea is to increase the par- 
ticle number in the initially unpopulated mp — Q spin 
component via spin dynamics transfer out of the other 
components. For this we first prepare a partially con- 
densed mixture of the |— 1) and states. The resulting 
dynamics can be divided into two main successive steps 
which are illustrated in Fig.|21as a) and b). 

The first process is that spin dynamics populates the 
mp = Q state by converting m^? = ±1 condensate atoms 
into mp = atoms according to |+ 1) + |— 1) <-> |0) -I- |0) 
0, 0, Q, 0|. Due to its density dependence spin 
dynamics is practically restricted to the condensed frac- 
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mp = —l mp = mp = -\-l 





FIG. 2: Scheme of the dynamics, (a) Spin dynamics transfers 
population to the mp = state which thermahzes almost im- 
mediately, (b) When all thermal clouds are equally populated 
and thus the critical particle number in mp —0 is reached a 
condensate arises and 'free' spin dynamics can take place. 



tions, resulting in the production of ?7if = 'condensate' 
atoms, which however immediately thermahze into the 
mp = thermal cloud due to collisions with all thermal 
clouds (Fig. We want to emphasize at this point 

that thermalization is the fastest timescale 50 ms) of 
our system and therefore spin dynamics (ss 1 s) is only a 
means to produce the new component. The redistribu- 
tion of constant total energy among more thermal atoms 
during this process leads as a side-effect to a decrease 
of temperature T. This is similar to 'decoherence-driven 
cooling' of the JILA experiment iTj which in contrast to 
our system did not involve conversion between different 
condensate components. 

As soon as the critical particle number in the mp = 
thermal cloud is reached the phase transition in the 
mp = Q component takes place and a condensate frac- 
tion emerges (Fig. From this moment on the ther- 
mal clouds are and remain equally populated and provide 
a constant temperature reservoir of the system. There- 
fore 'free' spin dynamics may take place between the spin 
components of the condensate fractions, i.e. at constant 
total number of condensed atoms but still in touch with 
the reservoir of finite temperature. Thus spin dynam- 
ics mainly determines the final mp condensate fractions, 
which are not as a rule equally populated in contrast to 
the thermal clouds 

The experimental setup (for details see pro- 
duces BECs in an optical dipole trap which provides a 
spin-independent trapping potential. The trapping fre- 
quencies are 27r x 890 Hz vertically, 27r x 160 Hz hori- 
zontally and 27r x 20 Hz along the beam direction. Spin 
dynamics is suppressed during preparation of the initial 
spin state due to the high magnetic offset field of 25 G 
which is subsequently lowered to a value of 340 ± 20 mG 
to allow for spin dynamics. After a variable hold time 
of 0..30s the dipole trap is switched off and the released 
atoms are spatially separated by a Stern-Gerlach gradi- 
ent. Finally an absorption image is taken in order to 
determine BEG and thermal atom numbers by a simul- 



taneous fit of three parabolas and three Gaussians for the 
three mp components. 
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FIG. 3: Measured condensate and thermal atom numbers for 
the different spin states (marked as 'exp.') as function of 
different hold times. The lines represent solutions of the rate- 
equation model for two different sets of spin dynamics pa- 
rameters denoted as 'simi' and 'sim2' (see text for numbers). 
The moments when the critical particle number for mp = is 
reached in the simulations are marked by vertical lines. 

Fig. El shows the experimentally obtained BEG and 
thermal atom numbers versus the hold time compared to 
simulations of the rate-equation model which will be pre- 
sented later. We start with an initial mixture oimp —H 
both in BEG and thermal fractions. The preparation 
process leads to a remaining population of < 10 % in the 
mp = state. 

The experimental data demonstrate all of the previ- 
ously introduced dynamics only modified by loss pro- 
cesses. First a mp = thermal cloud arises and grows 
until the critical particle number is reached after 5. .10s. 
Note, this is the moment of equal populations of all ther- 
mal clouds. Subsequently a mj? = condensate fraction 
emerges. The data between 5. .10 s suggest that the ex- 
act moment of phase transition varies from shot to shot. 
Indeed this moment crucially depends on spin dynamics 
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as will be discussed later. Finally spin dynamics leads to 
a steady-state which decreases due to loss processes with 
an experimentally observed relative condensate distribu- 
tion of 40..45% mp =±1 and 10. .20% mF=0 

In the following we develop a simple rate-equation 
model which reproduces the main experimental obser- 
vations. We do not intend to give a detailed and thor- 
ough simulation of finite temperature BEC which would 
be quite involved and is subject of current theoretical 

activities [ii HO, m m 

Rather a basic model from a 
experimentalist's point of view is presented to stimulate 
a vivid discussion of finite temperature effects in multi- 
component BECs and introduce a number of single pro- 
cesses which yield the observed behavior. The model is 
based on a set of 7 variables , , , N^^ , ,N^,T 
where with X = — ,0, -I- denote the atom numbers 
of mp ~ — 1,0, +1 in the condensate fraction and 
the respective atom numbers in the thermal cloud. T is 
the system temperature and assumed to be equal for all 
components. The equations of motion read 



X 
X 

t,th ' 



X 
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and include the processes thermalization {N^fi^, Tth), 
spin dynamics {Nq^^sp)^ one-body losses {N^j^ 
body losses (iVj^g^) and evaporation {N^ 
single effects as well as an aditionally introduced phase- 
space redistribution will be discussed in the following. 

The thermalization rate ^th quantifies the coUisional 
transfer of condensed atoms into the thermal component 



X 

0,th 
X 



and is obtained via the 



where Nt = + N° 
relation ^thNt — ^thnt which takes into account the peak 
density of the thermal cloud fit to convert the density de- 
pendent rate 74,1 into7f,j = 7thw3(m/(27rfcBT))(3/2). The 
temperature- and spin-dependence of ^th is neglected. 
The system temperature T decreases as the conserved 
total energy is redistributed among more thermal atoms 
and given as 

Tth = -TjthNo , 



with No ^ Nq + A^o + ^0 ■ The used value -fth 



10~^*m'^/s leads to a thermalization rate ■jthNo of ~ 13 
1/s for A'o = 45000 which corresponds to our experiment. 

Spin dynamics is implemented by a simple coupling of 
the condensate atoms due to the relation | — 1) + | + 1) ^ 
|0) -|- 10) with two reaction rates "fspi and 'ysp2 for forward 
and backward reaction l23l: 



± 

0,sp 



o,sp = -2%pi« + 27.p2^o"^o+ 



One-body loss occurs with rate 71 independently of the 
spin state and equally in the BEC and thermal cloud. 
The value used is 71 — 0.011 1/s and corresponds to the 
measured 1 /e-lifetime of 90 s limited by background gas 
collisions 
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For three-body loss [31 we assume a spin-independent 
process, ignore possible changes in statistical factors due 
to multiple components and obtain 



X 



with C3 = 7/6cl and C2 = 15'^/^{UTT)-\mLj /h^)^/^ . 
The loss rate used is L = 5.8 x IG^''^ hi^/s 25]. 

The evaporation process due to finite trap depth ksTe 
is implemented by a particle loss of the thermal cloud 
connected with a decrease of the system temperature. 
The temperature dependence of the evaporation rate 7e 
is neglected and the loss reads 

Energy conservation leads to a change of temperature 

We use the Euler method to propagate the equations 
in discrete time-steps of duration At (e.g. {t + At) = 
N^{t) + At). After each simulation step a phase- 
space redistribution is carried out. This is important 
to introduce quantum statistics into the equations and 
can be regarded as spontaneous condensation if the crit- 
ical density is reached. The critical particle number is 
calculated as Nc — g3{l){kBT / {hui))^ and the following 
condition is checked for X = — , 0, + 



If (7V/^ > N,) 



N,f{t 
N^{t 



At) = N^{t) + {N^'{t)~N,) 
At) = N, . 



This re-condensation step is related to a temperature 
change obtained by total energy conservation as 



T{t + At) = T{t) 1 



No{t + At) - Nojt) 
Nt{t + At) 



The thermalization step and the phase-space redistribu- 
tion cancel out in the case of thermal equilibrium result- 
ing in steady condensate fractions and constant temper- 
ature. Nevertheless these steps are crucial to describe 
the occurrence of the new thermal components and new 
condensate fractions. As thermalization is the fastest 
timescale of the considered system a step-like description 
seems to be reasonable. 
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Our rate-equation model reproduces all experimentally 
observed thermal features even with a reasonable quan- 
titative accuracy as shown in Fig. O The initial con- 
densate populations were chosen as N-{0) = N+{Q) = 
45000, N^{0) = 7000 and the thermal atom numbers 
as N-{0) = N+{0) = 90000 and N^{0) = 12000 and 
T(0) = 288 nK. Evaporation parameters are 7e = 0.015 
1/s and Te — 500 nK. Simulations for two sets of spin 
rates ^spi = 1-6 x 10"^ 1/s, ^sp2 0.4 x 10"^ 1/s (siml) 
and 7,pi = 2.4 x 10"^ 1/s, 7^p2 = 0.6 x 10"^ 1/s (sim2) 
have been carried out. Although these two sets of pa- 
rameters differ by only 33 % the resulting moment of con- 
densation varies by more than a factor of two (4 and 9 s 
respectively). Indeed we have to assume that there is 
a shot-to-shot variation of spin dynamics in our exper- 
iment as initial phases are not controlled. It has been 
theoretically shown [2^ that spin dynamics crucially de- 
pends on initial relative phases. Another influence on the 
spin dynamics rates may arise from shot-to-shot varying 
densities. 

In its simplicity the rate-equation model allows to ob- 
tain a clear physical picture of the dominant thermo- 
dynamical aspects but it lacks coherent spin dynamics 
which has been reduced to simple rate-equations. This 
procedure seems to be suitable for the discussed regime 
and may be applied to further problems in this con- 
text. Nevertheless the detailed treatment of shot-to-shot 
variations, coherent dynamics, excitations and phase- 
fluctuations of condensates 0| requires an extended the- 
oretical description. 

Finally we want to point out that the complementary 
case has been studied for = 2 of ^^Rb |12j, where 
spin dynamics (wlOms) is faster than thermalization 
(w 50 ms) leading first to a steady distribution of con- 
densate spin components which afterwards melt. 

In conclusion we have reported the experimental re- 
alization of a new regime of Bose-Einstein condensation 
in multi-component spinor condensates at finite temper- 
ature. A rate-equation model has been presented to dis- 
cuss the main thermodynamical aspects. The physics 
introduced here paves the way towards general new as- 
pects in multi-component quantum gas thermodynamics 
at finite temperature. 

We acknowledge support from the Deutsche 
Forschungsgemeinschaft in the SPP 1116. 
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